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GASIFICATION  EFFECTS  IN  THE  HETEROGENEOUS  IGNITION 
OF  A CONDENSED  FUEL  BY  A HOT  GAS 

M.  Kindelan  and  A.  Lilian 
ABSTRACT 
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gasification  flow  are  taken  into  account  in  obtaining  an 
integral  equation,  involving  a single  parameter  A,  which 
is  solved  to  describe  the  evolution  of  surface  temperature 
with  time.  A runnaway  in  surface  temperature  is  found  to 
occur  at  a well  defined  ignition  time,  which  is  calculated 
as  a function  of  A.  For  values  of  A above  a critical  value 
no  ignition  occurs  because  the  cooling  effects  of  the 
gasification  flow  dominate  over  the  effects  of  chemical 
heat  release. 
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Introduction 


Different  theories  have  been  pr  ^osed  to  explain  the  ig- 
nition mechanism  of  a condensed  material.  These  theories  can 
be  classified  as  condensed-phase,  gas-phase  and  heterogeneous 
ignition  theories  according  to  the  location  of  the  reaction 
responsible  for  ignition.  Experimental  results  backing  each 
of  these  locations  have  been  presented  in  the  past. 


In  this  paper  we  develop  an  asymptotic  analysis  of  igni- 
tion of  a condensed  material  which  is  suddenly  exposed  to  a hot 
ambient  gas.  These  conditions  occur,  for  instance,  when  a shock 
wave  is  reflected  by  a condensed  material  located  at  the  end  of 
a shock  tube.  We  will  consider  the  case  in  which  an  exothermic 
heterogeneous  reaction  takes  place  at  the  surface  of  the  material 
generations  gaseous  products. 


Anderson  et  al  . proposed  the  heterogeneous  ignition  me- 
chanism trying  to  explain  experimental  evidence  showi. g the  fea- 
sibility of  igniting  a condensed  fuel  by  contact  with  a powerful 

oxidizer  (hypergol).  An  analytical  description  of  this  model  was 

2 

carried  out  by  Williams  , using  Laplace  transforms  to  obtain  an 
integral  equation  relating  the  surface  temperature  to  time.  Al- 
though the  formulation  of  the  problem  included  a flux  of  radiant 
energy  absorbed  at  the  surface,  only  the  case  of  zero  flux  was 
numerically  integrated.  For  large  adtivation  energies  the  nume- 
rical results  exhibited  "runaway  in  surface  temperature"  at  a 
finite  time  (ignition  time). 


Williams'  analysis  has  been 

effects  of  in-depth  absorption  of 
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Closed  form  expressions  for  the  ignition  time  as  function 

of  the  parametersentering  the  problem,  have  been  derived  by 
5 

Waldman  using  both  a local  similarity  hypothesis  and  an  adapta- 
tion of  a method  proposed  by  Rae. 

Grishin  and  Ignatenko6,  calculated  the  first  three  terms 
of  an  expansion  of  the  surface  temperature  in  powers  of  the 
square  root  of  time.  The  resulting  surface  temperature  history 
exhibits  an  inflection  point  at  a time  which  they  define  as  the 
ignition  time.  However  the  expansion  for  the  surface  temperature 
does  not  converge  for  times  close  to  the  time  of  thermal  runaway. 

7 

Yung  and  Chung  used  the  Laplace  method  for  large  activation 
energy  in  an  incorrect  way  to  obtain  an  ordinary  differential 
equation  describing  the  surface  temperature  evolution  with  time. 

By  numerically  integrating  this  equation,  the  authors  obtain 
approximate  values  of  the  ignition  time,  although  the  analysis 
is  not  asymptotically  correct. 

Q 

Lilian  and  Crespo  have  ana 
process,  under  radiant  and  shock 
of  high  activation  energy.  They 
for  the  surface  temperature  invo 
parameter;  the  solution  exhibite 
time . 

The  present  paper  generalizes  the  results  of  the  preceeding 
analyses  for  ignition  under  shock  tube  conditions  by  including 
the  effects  of  convection  that  were  neglected  in  all  of  them. 

It  is  shown  that  in  general  the  convective  effects  are  important 
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in  heterogeneous  shock-tube  ignitions,  and  therefore  they  should 
not  be  neglected,  without  justification. 

Upon  contact  with  the  hot  gas,  the  surface  temperature 
raises  instantaneously  to  an  intermediate  "jump  temperature" 
which  is  a function  of  the  initial  temperatures  in  the  solid 
and  gas  and  of  the  ratio  of  thermal  responsivities . In  absence 
of  chemical  reaction,  the  surface  temperature  remains  constant 
and  the  width  of  the  conduction  layer  increases  with  the  square 
root  of  time.  However,  in  the  presence  of  a surface  reaction  there 
are  two  effects  that  tend  to  modify  the  temperature  at  the  sur- 
face. First,  because  of  the  gasification  flow  the  hot  gases  are 
blown  away  from  the  surface  and  cold  condensed  material  is  brought 
towards  the  surface  thereby  decreasing  the  surface  temperature. 
Second,  the  energy  released  by  the  exothermic  reaction  has  the 
effect  of  raising  the  surface  temperature. 

Two  different  regimes  are  observed  depending  on  which  of 
these  two  competitive  effects  dominates.  If  the  effect  of  the 
convective  flow  dominates  the  surface  temperature  continuously 
decreases  while  the  condensed  material  gasifies.  The  surface 
regression  rate  is  maximum  at  the  time  of  contact  with  the  hot 
gas  , and  thereafter  decreases  with  time  approaching  zero  for 
long  times.  If  the  effect  of  heat  released  by  the  reaction 
dominates,  the  surface  temperature  increases  with  time  thereby 
accelerating  the  reaction  rate.  This  self accelerating  behaviour 
causes  thermal  runaway  of  the  surface  temperature  at  a finite 
ignition  time. 

In  the  following  analysis  we  derive  an  integral  equation 
with  a unique  parameter,  to  describe  the  time  evolution  of  the 
surface  temperature.  This  parameter  is  a measure  of  the  effects 


4.- 


of  convection  and  has  a critical  value  which  separates  the  two 
different  regimes. 

2 . F ormulation 

We  consider  a one  dimensional  model  with  the  solid  occupying 
the  half  space  x < 0 and  the  oxidizer-containing  gas  the  region 
x > 0. 

Figure  1 is  a schematic  representation  of  the  process  and 
shows  the  effects  accounted  for  in  this  study. 

It  is  assumed  that  an  Arrhenius  heterogeneous  reaction,  of 
arbitrary  order  with  respect  to  oxidizer  mass  fraction,  occurs 
at  the  interface  between  the  condensed  fuel  and  the  ambient  oxi- 
dizer. For  convenience  the  origin  is  fixed  at  the  surface  of 
the  regressing  material  and  the  mass  coordinate  4)  is  introduced 
as  the  new  space  coordinate  in  the  gas  phase.  We  will  use  the 
well  justified  assumptions  that  the  work  associated  with  viscous, 
pressure  and  external  forces  is  negligible.  We  use  Fick's  law 
to  calculate  the  diffusion  velocities  and  consider  the  specific 
heats  to  be  constant  in  each  phase.  It  is  also  assumed  that 
the  density  and  thermal  conductivities  in  the  condensed  phase  as 
well  as  their  product  in  the  gas  phase  are  constant.  Under  these 
conditions  the  equations  of  conservation  of  energy  in  the  solid 
and  gas  phases,  and  the  oxidizer  concentration  equation  are 
respectively 
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where  all  the  symbols  are  defined  in  the  nomenclature.  The 
initial  and  boundary  conditions  are 

T(  0 , x)  = T(t,  — ) = Tq  , Tg(0,  *)  = Tg ( t » »)  = Tgo  , 

T(t,  o)  = T (t,  o)  (4) 

o 

Y Co , *)  = Y ( t > -)  = Yw  (5) 

_ x I + p A llfi.  = pvfcc  -c)(T  -T.  )~q!  (6) 

3x  |s  ^g  g 3^  s L 6 8 1 J 

(7) 

where  the  surface  regression  rate  v is  given  by 

v = BYn  exp ( -E/RT  ) (8) 

s s 

and  m is  the  stoichiometric  coefficient  in  the  reaction 

+ m [o]  ->■  products  (9) 

At  time  zero  the  temperature  in  the  gas  phase  is  suddenly 
raised  to  a higher  value  TgQ.  For  early  times  the  temperature 
gradients  at  the  surface  are  very  large,  and  therefore  the  right 
hand  side  of  equation  (6)  is  negligible.  Under  these  conditions 
the  problem  is  identical  to  the  heating  of  an  inert  condensed 
materM  under  shock  tube  conditions,  whose  solution  is  well 
known  in  the  case  of  constant  values  of  density,  heat  capacity 


and  conductivity.  We  derive  in  an  appendix  the  solution  to  the 
inert  problem  for  arbitrary  dependence  of  the  thermal  properties 
on  temperature.  In  the  analysis  that  follows,  the  quantity  p X 

© o 

is  taken  as  constant,  but  its  value  is  chosen  in  such  a way  that 
the  correct  "jump  temperature"  (temperature  at  the  surface  with 
inert  heating)  derived  in  the  appendix  is  obtained  when  the  inert 
problem  is  solved  in  the  assumption  of  constant  thermal  properties. 

Introducing  the  nondimensional  variables 
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in  the  solid  phase 
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where  is  the  jump  temperature 

•1  = To  ♦ <V  - V iTr 
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and  v^  is  the  value  of  v for  Tg  = and  Yg  = 
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equations 


e(o,  s)  = e(x,  °°)  = -r  » e (o,  ?)  = 0 (x,  °°)  = 1 

o o 
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3 . Analysis  of  the  ignition  process 

Upon  contact  with  the  hot  gases,  the  surface  temperature 
will  be  instantaneously  raised  to  its  jump  value.  However,  the 
e .^thermic  reaction  at  the  surface  introduces  two  effects  that 
modify  the  jump  temperature.  On  the  one  hand,  the  gasification 
flow  brings  cold  material  towards  the  surface  and  pushes  the-  hot 
gases  -way,  thus  reducing  the  surface  temperature.  On  the  other 
hand  the  heat  released  by  the  exothermic  reaction  causes  an 


J 
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increase  in  surface  temperature.  Ignition  will  or  will  not 
occur  depending  on  which  effect  dominates. 

It  is  observed  that  for  large  values  of  the  nondimens ional 

activation  energy  E/RT  , small  increments  in  surface  temperature 
o * 

of  order  RT^/E  produce  an  increment  in  the  reaction  rate  by  a 
factor  e. 


The  corresponding  increments  of  0 and  8g  from  their  inert 


values 


Oj/ r = erf ( S/2/7 ) 


(21) 


are  small,  of  order  RTT/lElT  -T,)|,  so  that  it  is  convenient  to 

1 go  l 

introduce  the  nondimens ional  variables 


x = 6(e-ex)/r 


xg  = s(eg-eIg)/r 


(22) 


= 6(Ya#  - Y)/[(Yw  ♦ m)  ar  (i+r  )/T] 


(23) 


a = t/  r . 


z = 6e/r 


(24) 


where 


B = Edj  - Tq)/RT' 


(25) 


is  a iiOndimens ional  activation  energy.  The  nondimensional  time 
o and  fuel  mass  fraction  y have  been  defined  so  that  they  are 
of  order  unity  during  the  ignition  transient.  The  analysis  that 
follows  is  based  on  the  assumption  that  the  parameter  B is  large 
compared  to  unity  which  is  equivalent  to  requiring  that. the 


T 
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reaction  rate  at  the  jump  temperature  be  much  faster  than  the 

gasification  rate  at  the  initial  temperature.  Thus,  it  is 

necessary  to  assume  that  3 is  large  in  order  to  have  a physically 

realistic  problem  in  which  n.  non  reacting  condensed  material  is 

suddenly  brought  to  vigorous  reaction  upon  contact  with  a hot 

oxidizer-containing  gas.  The  ratio  of  thermal  responsivities  T 

is  usually  small,  but  it  may  become  of  order  unity  for  large 

values  of  pressure  or  TgQ,  since  it  is  proportional  to  the  square 

root  of  p X . Thus,  the  parameter  6/  T is  at  least  of  order  6. 
g g 


Introducing  the  nondimensional  variables  (22)  - (24)  and 
taking  the  limit  6 °°,  equations  (14)  - (?0|)  simplify  to 
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exp  (-  + Xs) 


X ( 0 , z)  = x(cr,“)  = y ( 0 , z)  = y(a  ,»)  = 0 


x(o,  0)  = Xg(o>  0) 


(26) 

(27) 


(28) 


(1  + r ) exp(xs) 


(29) 
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(30) 


The  parameter  P appearing  in  equation  (26)  describing 
the  temperature  in  the  condensed  phase,  equals  a r (1  + T )• 


mm 


The  equation  that  describes  the  gas-phase  temperature  is  identical 
to  Eq.  (26)  with  P replaced  by  ar(l+r  )/  r 

Only  the  gradient  of  the  inert  temperature  appears  in  the 
effect  of  convection,  given  by  the  right-hand  side  of  Eq.  (26), 
and  no  convective  effects  appear  in  Eq.(27)  describing  the  mass 
fraction.  It  should  also  be  observed  that  in  the  limit  6 ■*  00 
the  equations  and  boundary  conditions  describing  the  temperature 
distributions  are  uncoupled  from  those  describing  the  mass  fraction 
Green  functions  may  be  used  to  write  the  solution  to  equations 
(26),  (28)  and  (29)  as 


x(°  > z)  = 


p o exp(x'  ) - 

- — / / exp(-z ' 2/4o ’ )G( z ,a ;z ' ,o  ' )dz  ’ da'  + 
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+ / ( 
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)do 
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where 


G(z,a;z,,a')  = exp[]-(z-z,)^/4(o-a,)]+exp£-(z+z,)^/4(o-a,)3(32) 

Evaluating  Eq.(3i)  at  the  surface  and  using  equation  (29) 

an  integral  equation  defining  x_  obtained,  viz., 

s 

o exp(x')  A a 

y = f da*  - ! exp(xi)da'  (33) 

8 o /it(o-o"  ) o s 

This  equation  describes  the  evolution  with  time  of  the 
surface  temperature  as  a function  of  the  single  parameter  A=a(r+r) 
which  is  a measure  of  the  cooling  effect  of  the  gasification 
flow.  The  range  of  values  of  A is  large;  when  the  overall  reac- 
tion responsible  for  ignition  has  a small  exothermicity  Q,  A 
may  become  of  order  unity. 
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' 

Equations  (27),  (28),  (30)  nay  be  used  to  obtain  through 
an  analogous  procedure,  an  integral  equation  relating  the  sur- 
face mass  fraction  of  oxidizer  to  time,  viz., 

o exp l x ' ) 

y = / - do*  (34) 

o /n(o-o ' ) 


The  first  term  in  the  right  hand  side  of  Eq.(33)  accounts 

for  the  heat  released  by  the  exothermic  surface  reaction,  and 

it  is  the  only  one  appearing  when  the  convective  effects  are 
4 

neglected  . However,  to  neglect  the  convective  effects  is  appr£ 
priate  only  when  A is  small;  that  is  when  the  nonaimensional 
heat  of  reaction  1/a  is  large  compared  to  unity  or  when  both 
the  ratio  of  heat  capacities  r and  the  ratio  of  thermal  respon- 
sivities  p are  small  compared  to  one. 


Since  in  Eq.(33),  the  term  associated  with  the  effect  of 

heat  released  by  the  reaction  tends  to  increase  the  surface 

temperature,  while  the  term  associated  with  the  effect  of  the 

gasification  flow  tends  to  decrease  it,  different  regimes  are 

observed  depending  on  which  of  these  two  competitive  effects 

dominates.  There  is  a critical  value  A =2  such  that  for  A > 2 

c 

the  effect  of  the  gasification  flow  dominates  and  therefore  the 
surface  temperature  continuously  decreases  while  the  hot  gases 
are  blown  away  from  the  surface.  For  A < 2 the  effect  of  heat 
released  by  the  reaction  dominates,  leading  to  a selfacceleratory 
reaction  rate  that  causes  runaway  in  surface  tempereture  at  a 
finite  ignition  time.  For  A = 2 the  surface  temperature  remains 
constant  at  its  jump  value,  and  the  temperature  distribution  is 
given  by 


erf  c C ) 

2 /o’ 


X 


Pz 

2 


(35) 
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When  A > 2 the  problem  is  similar  to  that  of  the  gasifi- 
cation of  a condensed  material  which  reacts  endotermically  at 
the  surface  with  the  ambient  oxidizer.  This  problem  has  been 
analyzed  in  Ref. (11)  where  an  integral  equation  similar  to  (33) 
is  derived.  For  very  long  times  such  that  x becomes  of  order 
6,  the  complete  convective  effects  have  to  be  retained  and  the 
Arrhenius  exponent  can  not  be  linearized.  This  regime  is  also 
analyzed  in  reference  (li)  where  a similarity  solution  is  derived, 
thus  completing  the  analysis  of  the  gasification  history.  When 
A < 2 , a finite  difference  version  of  Eq.(33)  has  been  numeri- 
cally integrated  through  a procedure  paralleling  that  used  in 
reference  9.  The  results  of  such  integrations  are  shown  in 
figure  (2)  for  several  values  of  A.  In  figure  (3),  the  nondi- 
mensional  ignition  time  a.  is  plotted  as  a function  of  A. 

When  A approaches  its  critical  value,  the  ignition  time  tends  to 
infinity . 

For  small  times  equation  (33)  yields 


X = /oT*  (2-A)tcr{l-A/2  + 2A(A-2)/3u}+.  . . 
s 


(36) 


This  expansion  has  been  previously  derived  for  the  case  A = 0 

Equation  (34)  has  also  been  numerically  integrated  to 
produce  the  time  evolution  of  the  oxidizer  mass  fraction  at 
the  surface.  These  results  are  shown  in  figure  (2). 


Discussion 

We  have  carried  out  an  asymptotic  analysis  of  the  hetero- 
geneous ignition  of  a condensed  material,  including  the  cooling 
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. effect  of  convection.  The  lowest  order  solution  for  the  surface 
temperature  deviation  from  its  inert  value,  in  an  expansion  in 
powerte  of  the  nondimens ional  activation  energy  6,  is  obtained 
by  solving  an  integral  equation.  A single  parameter  appears 
in  this  equation  to  measure  the  effects  of  convection,  and  when 
these  effects  are  dominant  it  is  found  that  the  gasification  flow 
moves  the  region  of  hot  gases  away  from  the  surface  and  no  igni- 
tion occurs,.  There  is  a critical  value,  A = 2,  of  the  convec- 

c 

tive  parameter  A that  separates  conditions  under  which  ignition 
occurs,  from  conditions  in  which  the  surface  temperature  conti- 
nucus'ly  decreases  with  time.  The  existence  of  these  critical 
conditions  has  not  been  previously  suggested,  since  this  criti- 
callity  is  associated  with  convective  effects  and  previous 
analyses  have  neglected  these  effects. 


Fob  values  of  A<2,  there  is  a thermal  runaway  at  an 

ignition  time.  t.  =t  o.  (A):  the  function  a.  (A)  is  repre- 
& ign  c ign  ign  * 

sented  in  Fig. 3,  and 


t 

c 


R2Tj(l+r)2  X c exp(2E/RT1) 
Q2  Y2n  B2  E2  p 


*37) 


where  T^  is  the  jump  temperature 

characteristic  time  is  inversely 

oxidizer  mass  fraction  raised  to 

the  heterogeneous  reaction.  The 

temperature  is  mainly  exhibited 

in  the  definition  of  t . Figure 

c 

dependence  of  the  ignition  time 


defined  in  Eq.(13).  The 
proportional  to  the  ambient 
twice  the  reaction  order  of 
dependence  on  the  initial  gas 
through  the  term  eXpy ( 2E/RT^  ) 
(4)  shows  e-sfclicf  telly  the 
on  the  initial  gas  temperature 


for  all  other  parameters  fixed.  It  is  observed  that  the  ignition 
time  decreases  rapidly  with  increasing  initial  gas  temperature. 
However,  when  this  initial  temperature  is  high  enough  as  to  make 


the  parameter  A close  to  2,  it  is  seen  from  figure  (3)  that  the 

nondimensional  ignition  time  cu  approaches  infinity.  In  figure 

ig  n 

00  it  is  observed  that  for  high  values  of  T , t.  increases 

* O 

with  the  initial  gas  tempa^ature  approaching  infinity  when  A goes 
to  2 . 

Our  analysis  also  shows  that  for  large  activation  energy, 
the  change  in  oxidizer  mass  fraction  from  its  initial  value  is 
small  and  therefore  it  can  be  neglected  when  analyzing  the  igni- 
tion process. 

This  study  is  applicable  to  the  case  in  which  an  endother- 
mic surface  reaction  produces  a gaseous  fuel  which  reacts  with 
the  ambient  oxidizer  immediately  after  gasification.  In  this 
case  the  heat  of  reaction  Q should  be  replaced  by  the  difference 
between  the  heat  of  reaction  andthe  heat  of  gasification,  as 
discussed  in  reference  (11) 
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Appendix 


The  ignition  time  of  a condensed  fuel  under  shock  tube 
conditions  has  been  found  to  depend  strongly  on  the  "jump" 
temperature.  Therefore,  it  is  important  to  compute  correctly 
its  value,  taking  into  account  the  dependence  of  the  thermal 
properties  on  temperature. 

The  jump  temperature  is  the  surface  temperature  instan- 
taneously established  after  contact  of  a condensed  material 

1 0 

with  a hot  ambient  gas.  Waldman  analyzed  this  problem  and 
derived,  t^e  correct  "jump"  temperature  by  assuming  the  Chapman- 
Rubesin  approximation  . In  the  following  analysis,  we  generalize 
his  results  by  allowing  an  arbitrary  temperature  dependence  of 
the  thermal  properties  in  the  gas  phase. 


The  equations  describing  the  process  are  equation(l), 

(2),  (4),  (6),  (8)  with  v = 0.  These  equations  describe  the  tem- 
perature distribution  in  the  absence  ff  surface  reaction  and 
gasification.  No  characteristic  length  or  time  appears  in  these 
equations;  therefore  it  is  easy  to  show  by  dimensional  analysis 
arguments  that  the  aquations  admit  a similarity  solution,  with 
the  surface  or  jump  temperature  T^  remaining  constant  for  t<0. 
Therefore  the  inert  solution  in  the  solid  is  still  given  by  Eq(21) 
if  its  density  specific  heats  and  conductivity  are  constant;  in 
the  gas-phase  the  energy  equation,  when  written  en  terms  of 

9 = T /T.  and  the  similarity  variable 
g 1 


Cp  A ). 

g g • 


CA1) 


takes  the  form 
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-L  J 7Ps\  li'Un  ^ = 

3n  | (Pgxg)1  an  J n 3n 


( A2  ) 


where  the  term  - 4t(30/3t)  has  been  left  out  of  the  equation, 
because  0 is  only  a function  of  n*  The  surface  boundary  condi- 
tions are  0(0)  = 1 and 


ii  = (1  . 

^ n /T  T1 

n=o 


(p  X c ) 
g.  g g 1 


(A3) 


obtained  from  Eq(  6 )•  At  infinity  0(°°)  = Tg0/T^. 

To  calculate  T^  and  the  solution  of  Eq(A2)  for  a given 

dependence  of  p X^  on  temperature  we  proceed  as  follows.  We  can 

integrate  Eq(A2)  to  calculate  0(°°)  = Tg0/T^  in  terms  of  (30/3n)ri  = 0 

or  equivalently  T^  if  Eq  (A3)  is  taken  into  account.  The  relation 

T = T (T . ) enables  as  to  calculate  T„  for  given  T 
go  go  1 i b go 

For  the  particular  case  in  which  p X_  is  assumed  to  be 
e g g 

proportional  to  the  temperature  raised  to  an  arbitrary  power  w. 


p X T 

g-g =(_&)« 


(A4) 


equation  (A2)  becomes 


2n  fi  + ±-  C0W  ii  ) = 0 

3n  3n  an 


(A5  ) 


with  boundary  conditions  0(o)  = 1,  3 0 / 3 n | fixed  . A Runge-Kutta 


' • 
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1 


1 


integration  scheme  has  been  used  to  solve  equation  (A5)  for 

different  values  of  w.  In  figure  (5)  the  ratio  T /T,  is 

go  1 


ed  as 

a funct 

ion  of 

30/3n  • To  us 

given 

conditi 

ons , g 

ive  values 

to  T 

e (5) 

T as 

a f unc 

tion  of  T, . 

One 

go 

1 

with 

the  giv 

en  val 

ues  of  T 

go 

to  d 

orresponds  to 

P 1 

g g 

independent 

of 

case 

the  clas 

sical 

result  given  in 

ered . 

In 

general 

the  p 

arameter  T 

shou 

Tj  - T0 

~T  - T. 
go  1 


(A6) 


where  is  the  correct  jump  temperature  derived  by  the 
method  outlined  in  the  previous  paragraph. 


These  results  may  be  easily  generalized  to  include 
an  arbitrary  dependence  of  the  heat  capacities  on  temperature 
as  well  as  for  variable  properties  of  the  solid* 


p-'i 
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NOMENCLATURE 


Pre-  exponential  factor 
Specific  heat  at  constant  pressure 
Binary  diffusion  coefficient 
Activation  energy 
Lewis  number 

Stoichiometric  coeficient 

Reaction  order  with  respect  to  oxidizer  mass  fraction 
Nondimensional  parameter  defined  by  aT(l+r)  in  the  con- 
densed phase  and  by  ar(l+F)/r  in  the  gas  phase. 

Heat  evolved  in  surface  reaction  per  unit  mass  of  fuel 
Universal  gas  constant 

Ratio  specific  heats  c /c 

& 

Temperature 

Time 

Surface  regressxon'.ca.te 

Exponent  of  the  dependence  of  p A on  temperature 

S g 

Space  coordinate 

Oxidizer  mass  fraction  : ? ’ i.  ; ; . ' i C ) 

Reduced  mass  fraction  defined  in  Eq.  (23) 

Reduced  space  coordihate  defined  in  Eq.  (24) 

Dimensionless  inverse  heat  of  reaction  c(T  -T. )/Q 

go  1 

Dimensionless  activation  energy  defined  in  Eq.(25) 

Ratio  of  thermal  responsivities  /(p  Ac ) / ( p A c ) 

S 

Dimensionless  parameter  a(f+r) 

Similarity  variable  defined  in  Eq.(Al) 

Dimensionless  temperature  (T-T^)/(T  -T^) 

Thermal  conductivity 

Dimensionless  space  coordinate  defined  in  Eq.  (11) 
Density 

Reduced  temperature  defined  in  Eq . (24) 

Dimensionless  time  defined  in  Eq . (12) 

Reduced  temperature  defined  in  Eq.(22) 

^ rx 

Mass  coordinate  p dx 

0 g 
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Ignition  time  t.  as  a function  of  the  initial 
ign 

gas-phase  temperature  for  particular  values  of  the 
other  parameters  of  the  problem. 
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Fig. 5-Ratio  of  the  initial  gas-phase  temperature  to  the  jump 

temperature  as  a function  of  the  nond imens ional  heat  loss 
to  the  condensed  phase  for  particular  values  of  the 

exponent  w of  the  temperature  dependence  of  p X . 


